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ABSTRACT 


An interactive package of computer programs has been developed for 
the analysis of time series data. The package, called the Time Series 
Editor, is designed around the Box-Jenkins' statistical methodology of 
time series analysis. The Time Series Editor was developed for time- 
shared use on the Controlled Program/Cambridge Monitor System (CP/CMS) 
but could be easily modified to accommodate other time-sharing systems. 
The Time Series Editor assists in data preparation, entry, analysis and 
diagnostic testing. Utilization of the package requires only a limited 
knowledge of the computer system with all required user responses 


prompted by the Editor. 





TABLE OF CONTENTS 


INTRODUCTION 


FPO CUP Ph O METHODOLOGY = = = = =- — — <- 


PROPERTIES OF STATIONARY PROCESS = = 


AUTOREGRESSIVE-MOVING AVERAGE MODELS 


HOMOGENEOUS NONSTATIONARY SERIES - - 


SEASONAL TINE SERIES SS = = = = 


ESTIMATION => = = = = = = -‏ ع م فير 


ENEE CHECKING = < < = — = = = 


FORECASTING - - - - - - - - - - - - 


THE EXECUTIVE PROGRAM = = - - 
DATA ENTRY AND PROGRAM OUTPUT - - - 
l. Offline Card Input = = = = = - - 
2. Numeric Keyboard Input = - - - - 
u EMSS SE Feyboard Input = > 
PLOT PROGRAM - - - ے >= ہے ے‎ => => - 


DIFF PROGRAM ہ‎ - - - => => - => - => - - 


AUTO PROGRAM = - - - => => => => > መ= ےہ‎ - 


ESTIMATE PROGRAM = = = = = = = = >= سے‎ 


FORECAST PROGRAM = መ = = = = = = = 


SIMULATE AND GENERATE PROGRAMS ፦ - - 


A. 


G. 


DESCRIPTION OF TRE TIME SERIES EDITOR 


A. 


BE 


H. 


EXAMPLE TIME SERIES ANALYSIS 


SUMMARY AND RECOMMENDATIONS 


L 


BEI. 


IV, 






APPENDIX A: User's Guide to Time Series Editor = - = = = - - - - 
APPENDIX B: Sample User Session and Program Output - - - - - - - 
APPENDIX C: Program Listings - - - - - - - - - => - - - - - - - - 
LIST OF REFERENCES = - - - - - - - - - - ہے‎ - - - - - - => => => => - - 
INITIAL DISTRIBUTION LIST - - - e ےہ‎ - 





D 
| 





LIST OF FIGURES 


Box-Jenkins forecasting method - - - - 
Card deck for OFFLINE READ - - - - - - 
Plots of series UU - - - - - - - - - - 
Plots of autocorrelations of series C 


Plots of partial autocorrelations 
of seriessC - - - => - - - - - - - - - 


Plots of autocorrelations of lst 
differenced series C = = - ےہ >= >= ے‎ - - 


Plots of partial autocorrelations 
of lst differenced series C ----- 


Plots of the forecasted values adjoined 
to the original series C with 80% 
confidence limits about the forecasts 


Plots of simulated series with 
model (1,1,0) of series C = = = => سات‎ 


6۰ 


Ue 


Figure 
Figure 
Figure 
Figure 


Figure 


Figure 


Figure 


Figure 


Figure 









LIST OE TABLES 
Auto and partial autocorrelations 
一 一 
Values of lst differenced series C 


Auto and partial autocorrelations 


of Ist differenced series C - - - 


Values of simulations for model 
(1,1,0) of series C ------- 


Table I. 


Table II. 


Table III. 


Table IV. 






de INTRODUCTION 


operations researchers, statisticians, economists, market- 
ing personnel, managers and many others are frequently faced 
with the need to analyze time series data. In most cases 
their objective is to discover patterns or recognizable be- 
havior in historical data that can be used to construct mathe- 
matical models of the time series from which forecasts of 
future behavior can be obtained. The importance of being able 
to forecast future behavior accurately cannot be overempha- 
sized. Whether the data be budget expenditures, populations, 
natural resource consumption, prices, demands, economic indi- 
cators, stock-market prices, manpower levels or whatever, 
decision makers concerned with planning for the future must 
base their decisions on their best predictions about the 
future behavior of the time series, 

Until the late 1960's, the analysis techniques were pri- 
marily those of spectral analysis with heavy application of 
harmonie analysis and mathematical transform theory. Because 
of the mathematical sophistication required by the spectral 
approach che analysis capability resided fairly exclusively 
in the hands of mathematicians and engineers. Consequently, 
many naive forecasting methods such as moving averages, ex- 
ponential smoothing and decomposition analyses were adopted 
by the majority of the decision makers. Since the late 1960's 


the statistical analysis of time series, embodied primarily 





in the methodology developed by Box and Jenkins [Ref. 2], 
has received widespread acceptance. Because the Box-Jenkins 
mF 15 descriped in a vocabulary more familiar to opera- 
tions researchers, statisticians, economists and managers, 
more and more business and government decision makers are 
building models from past data to use for planning into the 
future. 

Many algorithms and computer programs for performing the 
analyses required by the Box-Jenkins approach have been de- 
veloped and are readily available from many sources. One of 
the best sources is the collection of FORTRAN computer sub- 
routines which resides in the International Mathematical and 
Beecttistcical Library (IMSL) [Ref. 4]. The major problem with 
using the available computer resources lies not with any de- 
mctency Of the algorithms or the programs, but with the very 
nature of the Box-Jenkins approach. The Box-Jenkins method 
is an iterative approach which is described in Figure l. 


[See Wheelwright and Makridakis, ref. 7.] 


Postulate a general 

class of models 
Identify the model 
which can be tenta- Stage l 












Estimate parameters 
for the model tenta- 


+ | Stage 2 
l Diagnostic check: | 

MT No TS The model ۱ 
۱ adegua tea 
=- -=-= =-=- ہے ےہ‎ — = yt -Yes- - - - - = 


Use model to Sea 
enerate forecasts 3 


51 ۵۶ ه‎ 1l. Box-Jenkins forecasting method. 
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Figure 1 shows that the Box-Jenkins method is a multi- 
stage, iterative process. It begins with the postulation of 
general class of models which has been found, experimentally 
to be extremely rich. Thereafter, the modeling procedure 
continues as a trial-and-error process with several decision 
points where the analyst is required to select the next 
direction based on the best information available to him. 
Each stage of the process outlined in Figure 1 may consist 
Of Several steps, and, even with the existing computer soft- 
ware resources, the modeling process is usually very time 
consuming. For example, a typical Box-Jenkins time series 
analyst using a batch processing computer system with access 
to the IMSL library of subroutines might perform the follow- 
ing sequence of tasks: 

1. Prepare time series data. 


2. Plot and visually examine the time series looking for 
memstarionarity, trends, deterministic patterns, etc. 


۰۳۲ ۳ ۲ à program to call the IMSL subroutine that 
calculates the mean, the variance, the autocorrela- 
tions and the partial autocorrelations. 


4. Plot the autocorrelations and partial autocorrelations. 
This provides the major information needed for 
identification of the time series. 


Rise a program to call the IMSL subroutine which 
transforms the time series to adjust for seasonal 
patterns, nonstationary behavior or other behavior 
which deviates from that assumed by the class of 
models postulated. 


6. Repeat steps 2) through 4) using the transformed 
data. 


7. Review the statistical properties of the autocorre- 


fa ۵ ( partial aútocorrelations for tentative 
identification of the model. 


dech 






٣۰۰۰۰۰۶٠۰ ٣٠٠٠٠ ده ود‎ to call the IMSL program that esti- 
mates the model parameters and computes the 
residuals. 


ece a program to perform goodness of fit tests. 


10. Analyze the residuals following steps 1) through 9) 
just as was done with the original time series. 


ll. Refine the model using the information obtained by 
the analysis of the residuals. 


12. Repeat steps 8) through 10). 

13. When an adequate model is obtained, write a program 
to call the IMSL subroutine which forecasts and 
determines confidence intervals for future values 
of the time series. 

Between each pair of steps the user must manually intervene 
and make a subjective decision based on the information 
available. Thus, a great deal of user interaction is required 
in order to determine a mathematical model and a forecast 
denat iion. Even with rapid computer turnaround time, the 
process can easily consume a day or more of calendar time. 
This report describes an effort to alleviate some of tne 
problems involved with modeling time series using the Box- 
Jenkins approach. An interactive computer package which pro- 
vides easy user access to the computational computer sub- 
routines available in IMSL and similar subroutine libraries 
was developed. The package, called the Time Series Editor, 
was written for time-shared use on the Naval Postgraduate 
School's Controlled Program/Cambridge Monitor System (CP/ 
CMS). Since all programs except the executive routine are 
written in FORTRAN, the Time Series Editor could be easily 


modified to accommodate other time-sharing systems. The 


Time Series Editor assists the user in data preparation, 
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Eu Oe construction, diagnostic testing and forecast- 
M re package requires only a limited 
knowledge of CP/CMS. In fact, with the User's Guide pro- 
vided in this report, a complete Box-Jenkins time series 
analysis can be performed in a short time by even a naive 
spmputer user. For this reason, the Time Series Editor 
should be valuable as an instructional aid for laboratory 
use in a time series class. 

A brief description of the Box-Jenkins methodology is 
given in Chapter 2 to serve as a point of reference for the 
remaining material. Chapter 3 contains descriptions of 
each of the programs contained in the Time Series Editor. 
The use of the Time Series Editor is illustrated with an 
example time series which is given in Chapter 4. 

Chapter 5 contains a summary and recommendations for 
additions to the Time Series Editor. A User's Guide which 
includes an explanation of CP/CMS sufficient for utiliza- 
tion of the Time Series Editor is given in Appendix A. 
Sample user sessions, sample outputs and complete computer 


listings are also included in appendices. 
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L PS BORZIENKINS METHODOLOGY 


In this chapter a brief description of the Box-Jenkins 
time series modeling methodology is given. For a more de- 
tailed discussion the reader is referred to the texts by 
Anderson [Ref. 1], Box and Jenkins (Ref. 2], Pindyck and 
Rubinfeld [Ref. 6], and Nelson [Ref. 5]. The material pre- 
sented here is included primarily to serve as a point of 
reference for the program descriptions that follow in later 
chapters. It is also included here to aid the user in 
understanding the computer output and the questions asked 


by the Time Series Editor. 


A. PROPERTIES OF STATIONARY PROCESSES 

A discrete stochastic time series is a set of observa- 
tions ሃ1፣። با ۰۰۶+ مولا‎ generated sequentially in time by a 
SST jointly distributed random variables; i.e., the data 
۷ ۰۰۰۶ You represents a particular realization of a joint 
prepapmimey distribution f (Y 1 Y oro." Ym)- A future obser- 


yation, can be thought of as being generated by a condi- 


YT+k 


7 گگ probability distribution function He‏ ا 3 دم 
given the realization through time T. The stochastic pro-‏ 
cess which generates the time series is said to be station-‏ 
ary if its properties are unaffected by a change of time‏ 
Sona nn i S it the joint probability distribution asso-‏ 


ciated with m observations ل‎ Vi made at any 


1 2 m 
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Sen of times ty Eure... p. 1S The same as that distribution 


associated with m observations y made 


E +78508 Yt tk 


cu times رر مت‎ t tx. 


mera precess 15 Stationary, the probability distribution 
£(y,) is the same for all times t. Thus the process has a 


constant mean 


CO 


u = E ly, = J senay 


- OO 
which defines the level about which it fluctuates, and a con- 
stant variance 


c? - EL (y =H) È) di (y-w) ^£(y)dy 


— CO 


which measures its variability about the mean level. Since 
Ene probability distribution E (yo) is the same for all times 
t, the mean and the variance can be estimated by the averages 
taken over time: 


N 
ye = 


2 | 


à 2 
and 5° = ES (y,-y) 


t=1 


Ale 


ሃኑ 
t=1 


The stationarity assumption also implies that the bi- 


39 ط۶1 3636 1ط‎ 7 toys Iove ) is the same for all times 
1 2 


[ana ch that lt = ISSO me tant The autocorrelation 


1 2 2 1 
at. lag x, Pur is defined as: 


0 


Et ty u) (Y 


0 
k o? 


ር+ጄ - 


This is estimated by the time average: 


H 
Il 
=] 
I 
PN 


= = 2 
UD UY ٤ 


ct 
ዘ 
a 


1882 








3۰۰۰٠٣٠٠٠٠۶٠۷٠٦٠ antocorrelation function vs. the lag k, 
ር... the correlogram, is very useful for the purpose of 
۶۳۲ 2۳11 ۲:9 Lf a process is stationary and for identifying 
the appropriate model. 

Another function which is important for purposes of 
identifying the appropriate linear time series model is the 


D. s 


^ 
a : : ۲ = $ + 
partial autocorrelation function Let Y, bot Pi ديعلا‎ 2 


where the b's are the least squares‏ ط000 


Yr+2 mS 


estimates of the linear regression coefficients (B's) in 


the model 


vee = + 517 +...+ ß 


+ 
o tl k-1Ytek-1 t 


Let ZL be the residual of I after removing the linear effect 


of Erom y Laes 


Ycekrl]'tt]Ye.k-1 c? 


^ 


2e — Ye. "Ze: 


The partial autocorrelation of lag k, denoted Pr ris defined 


to be the simple correlation of lag k for the adjusted series 


B, AUTOREGRESSIVE-MOVING AVERAGE MODELS 
The general class of models postulated by Box and Jenkins 
for stationary time series is the class of linear models 


defined by: 


16 


where ከ is a sequence of observations from a white noise 


2 
process (Ela, ] = 0, var la, ] A and E[a, a, |! acOCror all 


; s sess are p+q+1l parameters that 


are to be estimated from the data. The model above is 


E> 0) and 0 ٢۳٢ o 
called a mixed autoregressive-moving average (ARMA) model of 
order (p,q). If q=0, the model is called an autoregressive 
model of order p, d and if p=0, the model is called a mov- 
ing average model of order q, RA Thus, the general linear 
model of Box-Jenkins represents the current observation ee 
as a weighted sum of past observations and present and past 
random shock terms. The model is usually expressed in abbre- 


viated form using operator and transfer function notation: 


$(B)y, = 8, * O(B)a, (25 
where Q(B) = 1-0.B-0 Red. cd 5” O(B) » 1-8,B-...-8 B and 
1 2 p /” 1 q 
B is the operator defined by B تو‎ In order to guaran- 


É 
tee stationarity lt is necessary that the autoregressive 


parameters satisfy certain conditions. The conditions can 
be summarized by stating that all roots of the polynomial 
equation Q(B) = O (treating B as a dummy variable) must lie 
outside the unit circle. 

The tentative identification of the appropriate member of 
the general class (the identification of p and q) is accom- 
plished by comparing the sample autocorrelation and partial 
autocorrelation functions of the given time series with the 
2 (۰ ۰1 ۰ ۶۲۲ 9 0۲ ۲ 61 2 ۰10۳ 2۳ partial autocorrelation func- 


tions of members of the general linear class. For most 
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stationary time series an adequate fit can be found in a 


model with p and q relatively small, say three or less. 


er HOMOGENEOUS NONSTATIONARY SERIES 

In many cases the time series of interest is not station- 
ary. Instead, the probabilistic structure of the process 
which generates the time series may change with time. For 
example, there may be some sort of trend or seasonal pattern 
in the time series, If the process does, nevertheless, 
exhibit behavior which is somewhat homogeneous then the 
original time series can often be transformed into a station- 
ary series that can be described by an ARMA model. 

A time series is said to be homogeneous nonstationary of 


order d if w = Na 


5 is a stationary series. Here À denotes 


七 


the differencing operator: 


= (1-B)y, 


That a time series is nonstationary is indicated by a plot 

of the time series itself (e.g. a nonconstant mean) and by 
Mie autocorrelation function. Characteristic of the correlo- 
gram of a nonstationary series is the very slow damping out 
of the autocorrelation. When this property of the correlo- 
gram is observed, the user should difference the series one 
time and compute the correlogram for the series — Ay E 


me user snonla continue to difference the series until the 


resulting series appears stationary or until the procedure 
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appears not to improve the series. The transformed series 
(o Az, is then modeled as an ARMA model. The resulting 


model, in terms of the original time series is: 


d = 
d(B)V حر‎ D. + 9 ة (ظ)‎ ۲ (3) 


In this form of the model, the transfer function Q(B) is 
assumed to be stationary; i.e., all roots of Q(B) = O are 
murside the unit circle. It is sometimes written in the 


equivalent form 


ህ (5) ሃ.. = 7 + 9(B)a, (4) 


where W(B) = 6 (8) V2 - d (8) (1-8) © and, clearly, WB) is not 
a transfer function of a stationary series (W(B) = 0 has 4 
roots on the unit circle). The ARMA model for the differ- 
enced series is called an autoregressive integrated moving 
average (ARIMA) model of order (p,d,q). For the purpose of 
distinguishing between the two forms of the ARIMA model, 
equation (3) is referred to as the differenced form and (4) 


as the undifferenced form. 


EIS SEASONAL TIME SERIES 

Seasonality is defined as cyclical behavior that occurs 
on a regular calendar basis. For example, a highly seasonal 
time series would be the sales of Christmas ornaments which 
exhibit a strong peak every December. Rainfall, crop yields, 
livestock production, energy consumption, and many other 


time series that are influenced by the weather all exhibit 
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seasonal patterns. Seasonal patterns are often easy to spot 
simply by observing the time series directly. However, many 
times, if the variability in the time series is large, 
seasonal patterns will not be distinguishable from the other 
fluctuations. Recognition of seasonality is important since 
EN rovides information that can aid in modeling and fore- 
casting. The autocorrelation function makes recognition of 
Seasonal patterns easier. 

Suppose, for example, that a monthly series has an annual 
Seasonal pattern. Then the realizations should show some 
special correlation with other realizations which lead or lag 
by 12 months; i.e., there should be some correlation between 
and y These corre- 


y, and y 


t 1+12' Ye +24" Ze and Y.,3gr etc. 


lations should manifest themselves in the autocorrelation 
function which should show peaks at k = 12, 24, 36, etc. 

The Box-Jenkins modeling approach for seasonal (nonsta- 
tionary) time series is to first transform the seasonal series 
to a new series which is stationary. This can often be 
accomplished by taking seasonal differences acy, defined as 


follows: 


5 
À y, 3 y, - E را کے‎ 


rt 


n A ሀ) 


ው dal u ES 
ج‎ ሷ. (ል. AJ > (1B) Y. 


ct 


۱ : 4 ۱ 
The transformed time series ل ما‎ (W = A SY y) 1 5 then analyzed 


È 
as a stationary time series, (It may be necessary to perform 


more than one seasonal and/or differencing transformation 


Remo resulting seczies às Stationary.) Suppose that 
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the resulting ARMA model for the transformed series 


is 


b(B)w. = > - 9 )8( ره‎ 


Mae model for the original series is then 


Seel = 
$ (B) (1-B') Dea d ES 9 (5) 5. À 


E. PARAMETER ESTIMATION 
Suppose the series has been tentatively identified as 


an ARIMA (p,d,q) model: 


d 
$(B)A y, - 9, * 6(B)a, (5) 
ው = 5 p ም n E 7 8 
where Q(B) = 1 dB — $ 5B and O (B) H 9,B DEE 0 . There 
2 
are p+q+2 unknown parameters, ]9 ۰۰۰/۵۱9۵۶۱۰۰: چم‎ 


that must be estimated. The Box-Jenkins procedure separates 
the estimation problem into two parts. First estimates are 
obtained for the autoregressive-moving average parameters 
9 and 9, and the estimates are then made of 7 and 1 which 
are functions of the ARMA parameters. The usual procedure 
is to select those parameter values ó and 6 that minimize the 
sum of squared errors. Let GE Re and 

d 


E - A "ከ ከመም Then, it can be easily shown that و‎ E [w, ] 


and that (5) can be rewritten as: 
$ (B) (w Hg? = 0 )8( 3 


EE ut - 
or aL = 0 (B) $ (B) (w,-4,,) 
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Let w, the sample mean, be the estimate of EI due 


is usually 0) and let 


£ > 


- 8 ^ (8) 6 (8) (w,-w) | (6) 


人 ^ A 
Let 5(9,9) = La š 
15 
ሓ 人 
The objective is to select those parameters Q and 0 that mini- 
A A ۱ ۱ ۱ 
mize S(Q,0). Since the equation S is nonlinear in the para- 


meters, iterative search methods must usually be used in the 


- SE 
minimization. An estimate of Va 1S provided by 


8 " - S($ 8)/n-p-q 


H DIAGNOSTIC CHECKING 

After the model has been tentatively identified and para- 
meter estimates have been obtained, the next task is to test 
whether or not the original specification was correct and 
the model is adequate. The process of testing the model 
takes many forms, but usually involves at least the following 
two steps: 

l. Generate a simulated series from the estimated model 
and compare the Simulated series and its autocorrela- 
tion functions with the original series and its 
respective autocorrelation and partial autocorrela- 
tion functions. The comparison is primarily subjec- 


tive. 
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2. Compute the residuals of the estimated model from 
(6) and compare the properties of the residuals with 
those assumed for the shock terms of the actual pro- 
cess. The residuals should be normally distributed 
and uncorrelated with each other. There are many 
quantitative statistical tests that can be used to 
test hypotheses about normality and zero correlation. 
A plot of the autocorrelation and partial autocorrela- 
tion functions of the residuals provides not only a test of 
whether or not the residuals are uncorrelated, but, if they 
are correlated, the plots suggest modifications to the model. 
For example, suppose the model was tentatively specified as 
the ARMA(1,1) model: (1-0.5B) (y,-2) = (1+0.7B)a, and the 
autocorrelations and partial autocorrelations of the resid- 


uals suggest the model 


(1-0.3B)a, = u, 


where the u's are white noise (uncorrelated with variance 
2 eg 
e ). Then, the next model entertained for the original 


Eme series should be the ARMA(2,1) model: 


(1-0.3B) (1-0.5B) (y,-2) = (1+0.7B)u, ; 


Gr FORECASTING 

The objective of forecasting is to predict future values 
with as little error as possible. The criterion most often 
used for selection of the best forecast is that forecast 


which has minimum mean square forecast error. Thus, if 
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7+ represents the forecast at origin T of the value Zait? 


the objective is to select VR) so that 


A 2 
ን ከ 700)) | 


is minimized. This forecast is given by taking VE, as the 


conditional expectation of وب لا‎ * 


The forecasts can be easily generated recursively from the 


Mathematical model utilizing the fact that 


E [y .] = DIE for ge (T is current time) 


and 0 AE ET 


Ela, $ 


For example, suppose the estimated model is: 
2 
(1-0.5B + 0.6B IAR = (1+0.3B)a, 


Or, equivalently, 


፦ ۰ — ۰ + + = 
Ge 0 ہ لاڈ‎ 0.6y a 0 3a, _ 


EE è Í 
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Dagro and a = 0.2. The forecasts 


where Yioo ^ 1.4; 


Yoo 100 


o £ Yıoı' Y102' and لآ‎ made at time t = 100 are found as 


follows: 
Y1001D = Elvyga ligo! Vogt erva] 
= وم و۵ +وو 0۰6۷و و لا5 ۰ 0] ظ‎ ۲۵۰34 pg) 
= 9-3Y 100 - 0 ۰6۷ وم 0۰321 + وو‎ 
Y100!1) --ے‎ 6 
= 8 ]0 ۰5۷ ون و‎ 9-6y)ootayo2t0-3a,.1] 
= 0.5Y 70011) = 0.6y1oo 
Y100!2) = -0.76 
ለየ a ge 100" 
اا ل‎ 
= -0,48 
Let e, (2) ہج‎ Y (X) be the forecast error X periods 


(8) 


ahead. It can be shown that ep (2) is given by 


en (2) 5 Bing Qt Mame gua te tM gana 


where the weights Va are determined from 


Eee) ToB) 


The variance of the forecast error is given by 


2 
jo" 


2 2 
Ele, (2)] = (1 + um J a Voy 3 
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CO a confidence interval of z standard deviations 
around a forecast AL periods ahead would be given by 
2-1 


ር; (ሃሐ(5)) = Y,U) t z( €», 


on EO) 
J= 


272 
y 1 
J 

Note from expression (8) that the one-step ahead forecast 


error, eS, 0: is simply a l.e. 


T+] " 


TEI 


This explains the common use of the word residual to refer 
come random shock terms. Also, from expressions (9) and 
(10) it is clear that the forecast error variance is a non- 
decreasing function of the length of the forecast period ዲደ. 
Thus, the confidence bands must get wider as the forecast 


period gets larger. 
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LET. ۲۳ ۰۱۳ TON OE THE TINE SERIES EDITOR 


In this chapter, descriptions are given of each program 
in the Time Series Editor that interacts with the user. 
There are eight separate program modules that are currently 
١10 111946446 in the Editor. In addition, there are a few other 
programs that reside in the Editor, but which are completely 
ከ ሚግ 951 ۵1۳1 to the user, The latter programs serve useful 
Wwwéeledohns in conjunction with the other eight modules, but 
they are not described here. 

No attempt is made to describe the actual mathematical 
calculations or algorithms that are performed by the programs. 
Rather, the objective is to give general descriptions about 
what each program can do for the user and how the user in- 
teracts with the programs. The eight programs described in 
this chapter are TIMESER EXEC, PLOT, DIFF, AUTO, ESTIMATE, 


FORECAST, SIMULATE and GENERATE. 


A. THE EXECUTIVE PROGRAM 

The Time Series Editor contains a master program, called 
INES ER BAEC, that provides file control of all of the other 
programs, controls input and output, takes care of the neces- 
sary CP/CMS protocol and provides instructions to the user 
as to what is in the Time Series Editor and how each program 
can be used. TIMESER is written ina special CP/CMS Exec 


language. TE 15 Enhe only program in the Editor that ds not 
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in FORTRAN, and, consequently, it should be the only‏ ده كان 
program that would need modification if the Editor were to‏ 
be adapted to another time sharing system.‏ 

After the user has logged into CP/CMS and linked to the 
files containing the Time Series Editor (a user's guide for 
this is given in Appendix A) the entire Editor package is 
Made available to the user by his command, TIMESER EXEC.” 

On entry of this command, the EXEC provides a guided tour 
Enough the Editor. It tells the user what the Editor can 
provide; it asks the user what tasks he wants to do; and, on 
the basis of the user's answers, it instructs the user as to 
what data is required and how it must be entered. When the 
user selects an option for execution, the EXEC loads the 
appropriate program(s) and automatically manipulates any re- 


quired input and output files. 


B. DATA ENTRY AND PROGRAM OUTPUT 

Whenever data is required, the user is prompted by either 
the EXEC or the program module being executed. In most cases, 
the necessary user response is a single alphanumeric charac- 
ter input during execution by keyboard. However, in some 
cases, the amount of data required is too bulky for keyboard 
entry, and the data is entered more efficiently offline via 


cards or tapes. Similarly, most of the output is typed out 


EE MaS e coma executive routine, called TS EXEC, 
that can be used by the more experienced analysts who wish 
to suppress some of the user instructions. This abbreviated 
program provides the same basic services as TIMESER EXEC. 
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PAGINE ات‎ ۲ user's terminal, but, in some cases, the output 
1 9۳1 ۳ offline for conservation of time and to provide 
a hard copy of the results. 

Detailed descriptions are given of the input and output 
requirements of each program in the individual program write- 
ups. However, there are some general principles of data 
ክክ ያ. that apply for all programs. These are described in 
this section. 

۱۳ OfFline Card Input 

When a large volume of data, such as a time series 
of a hundred or more observations, is required, the data can 
be entered more efficiently through mechanisms other than 
keyboard entry. Keyboard entry would be not only much 
slower, but also more likely to contain errors than other in- 
put mediums. Thus, the Time Series Editor requires that the 
time series data be entered offline via cards. The data are 
read offline and stored in the user's file FTO2FOOl which is 
read automatically when required by the Time Series Editor. 
An example of an input data deck for use in the Editor is 


shown in Figure 2. 


29 





10th 2 ۹ ۹۹ٗ 4 07٤ 
| 1 











09 7706 2ی .27-0-7۰05 26268 


(5 observations per card, f 5E1546) 





250 
(number of observations) 


mol A | 


/ OFFLINE READ FILE FTO2F001 





aaaa is the user's identification number 
(the number starts in col. 13) 


Figure 2. Card deck for offline read, 


3:0 


ra 





2 Numeric Keyboard Input 

When the user is prompted to enter numerical values 
such as the number of observations, the number of parameters, 
parameter estimates or starting conditions, he must enter 
those values before program execution can continue. The user 
should enter the data according to standard FORTRAN practice. 
That is, integer data should be entered (without decimal) for 
counts and names beginning with the letters I through N; float- 
ing data should be entered with decimal point for all other 
variables. Because a typed decimal point overrides a float- 
۳ ۶۱۵۲ it is not necessary for the user to concern himself 
with the format for floating data (the user is never asked to 
enter more than a single observation on a line). However, 
some care must be exercised when entering integer data because 
integer data must be right justified in its format field. 
The user is told whenever the integer format is anything other 
Ea Il. For example, suppose the user wants to analyze a 
time series having only 65 observations. If the program that 
he is executing requires the length of the time series, the 
user will receive the following request: 

MIRER T CENGTH OE THE TIME SERIES, L, VIA I3. 
The user should then enter: 


Gol. 123 
b65 (b represents a blank space) 


If the blank were omitted, the program would read the length 


as 650 and many problems would occur. 
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a Alphabetic Keyboard Input 
The Time Series Editor often requests the user to 
Er spoma with alphabetic input. For example, it may ask him 
a question that requires a yes or no answer, or it may ask 
him what option he wants to execute. The editor has been 
programmed to read only the first letter of the user's re- 
sponse. Thus, he need only enter a single letter for each 
Such inguiry. For example, he should enter Y for yes, N for 
۰۰. cor the PLOT option, E for the ESTIMATE option, etc. 
4. Output 
Most of the results are written out right at the 
user's terminal. In some cases, however, the output is prin- 
ted offline to conserve time. Such results as plots and 
transformed time series are written onto various files 
END م۳ ول‎ 9 001, FTOIFOOl) and the EXEC program prints 
them offline under the user's identification number. The 
files can also be printed out at the user's terminal at his 


request. 


am PLOT PROGRAM 

The PLOT program plots any given time series which re- 
519555 in file FTO2F001. Other than the time series, which 
is entered offline, the program requires only that an identi- 
fication title for the plot be entered by the user during 
execution. 7۲۲ ۳1 ۳ ۶ automatically printed offline and is 
SS ound in file FTOSFOOL. The PLOT program uses the sub- 
routine PLOTP in the IBM Scientific Subroutine Package 


1۱ ۲ عه‎ (SSPLIB). 
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ES DIFF PROGRAM 

The DIFF program performs seasonal and/or nonseasonal 
differences and it allows the user to transform a given time 
series using logarithmic or exponential transformations. The 
program requires the given time series to reside in file 
۳۰ 02۳001, and it requires the user to input the following 
Wwéormatlon during program execution: 

(1) Type of time series (seasonal, nonseasonal) 

(2) Order of differencing (for stationarity) 

03) Length of seasonal period (for seasonal Series only) 


(4) Type of transformation (none, exponential, 
logarithmic) 


(271 Transform parameters (if log or exponential trans- 
form is desired) 


(6) Yes or no responses to questions about plotting of 
autocorrelations and partial autocorrelations of 
transformed/differenced series. 

The transformed/differenced series is written out into file 
۳۰۲ 0۵ 3۳5001 and can be printed out at the user's terminal if 

he requests a printout. Plots of the transformed/differenced 
series, its autocorrelations, and its partial autocorrela- 
tions can be printed offline. 

This program is used to transform a given series that may 

be seasonal and/or nonstationary or nonhomogeneous into a 
stationary series of the type possible for analysis using 


the Box-Jenkins methodology. The program utilizes the IMSL 


subroutine, FTRDIF. 


E. AUTO PROGRAM 
AUTO calculates summary statistics from a given time 
Series. The summary statistics include the sample mean, 
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Variance, autocorrelations, and partial autocorrelations of 
MESS one through 20. All of the statistics are printed out 
at the user's terminal, and, in addition, plots of the auto- 
correlations and partial autocorrelations are printed out 

St Eline. Furthermore, the plots reside in the user's file 
FTOSFOOl so that they may also be printed at the user's ter- 
minal. That is, however, a very time consuming process. 
AUTO provides the major information about stationarity, 
seasonality and model identification. With the Box-Jenkins 
procedure, the second moments (autocorrelations and partial 
autocorrelations) are the primary tools for tentative model 
۱۳۳۱۵ ۱۱۱ 1 ۳ 102۱ . 111 594492921 ,ذه‎ AUTO 15 used in the diagnostic 
checkout phase of model building to test the residuals. 

The 七 Ime series (öriginal, transformed or differenced, 
or residuals) must reside in file FTO2FOOl. (When a series 
is transformed or the residuals are calculated, those values 
are automatically stored in file FTO2FOO1 temporarily for 
further analysis.) The program uses the IMSL subroutine 


— 


Ri ESTIMATE PROGRAM 

After the user has tentatively identified a model, the 
Editor program ESTIMATE should be executed to calculate maxi- 
mum likelihood estimates of the model parameters. It esti- 
mates the autoregressive parameters, the moving average 
parameters, the constant term and the variance of the shock 
terms. It also determines the residuals which are so impor- 


eanesror che diagnostic checkout phase of model building. 
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AMATE USES the IMSL subroutine FTMAXL to estimate the 
parameters. 

ESTIMATE needs the time series in file FTO2FOOl prior to 
execution. During program execution the user will be prompted 
with the following requests: 

1) Enter number of autoregressive (AR) parameters. 

2) Enter number of moving average (MA) parameters. 

3) Enter the number of differences. 

ESTIMATE then provides the following output: 


1) Estimated AR parameters for the undifferenced form 
of the model. 


2) Estimated MA parameters. 

3 ) Estimated overall MA constant and white noise variance. 
4) Auto- and partial correlations of residuals. 

2 Pilots of (4). 


6) Values of residuals printed out offline and residuals 
stored in file FTO3FOOI. 


7) Chi-square goodness of fit value for estimated model. 


G. FORECAST PROGRAM 

The Editor program, FORECAST, uses the estimated mathe- 
matical model to compute forecasts of future values of the 
time series. It also computes (l-a) 100% probability limits 
for the forecasted values. The program utilizes the IMSL 
subroutine FTCAST. 

The time series must reside in file FTO2FOOl before exe- 
cution begins. During execution, the user is required to 


enter the following inputs: 
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2 ) Number of AR parameters 

3) Estimated values of AR parameters 

4) Number of MA parameters 

E) Estimated values of MA parameters 

6) Overall MA constant and white noise variance 

7) Maximum lead time for a forecast 

8) Order of differencing in model 

ED) Significance level for forecast confidence limits 
The program then provides the following output: 

1) AR parameters in undifferenced form 

for lead times L = 1,2,...,max‏ دا و 9 آلا 


3) Deviations from each forecast for the (l-a) 100% 
confidence limits 


4) Plots of forecasts and corresponding deviations joined 
with the original time series. This is plotted 
Offline- 

H. SIMULATE AND GENERATE PROGRAMS 

Because of their similarities the Editor programs SIMULATE 
and GENERATE are described together. GENERATE allows the user 
to generate a time series from any ARIMA model that he 
specifies. The user must identify the model and give values 
for its parameters and starting conditions. The program takes 
the specified model, generates random noise terms, and 


calculates as many values of the time series as desired. 


ee the identified model was ARIMA (1,1,0). The 
differenced form is (1-$,B) Ay, = Um +a. The undifferenced 


| E 
form is (1-(1+9,)B+0,B%)y, = D. ee የ ር ከ1 0796 


E 
(1-6, B) ۱5/۸ ۱ . 
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S program is useful for purposes of classroom instruction 
for the generation of a wide variety of time series examples 
men model identification. It is also useful for the diag- 
nostic phase of model checkout. A time series can be 
generated from the estimated model and its properties can be 
compared with those of the original series. If large dis- 
crepancies occur the estimated model may be inadequate. 

The program requires the user to enter the following in- 
pecs during execution: 

T) A randum number seed 

2) Number of AR parameters (undifferenced form) 

3) Number of MA parameters 

4) Length of time series 

5) White noise variance 

6) Values of AR and MA parameters 

7) Initial starting values 
The program output is the generated time series. 

The SIMULATE program provides the capability of generat- 
ing any number of simulated time series'. This is useful for 
predicting what might happen in the future and to demonstrate 
that, even within a given model, the actual observed time 
series' can differ substantially. This program uses GENERATE 
but also requires as input the number of simulated series the 
user wishes to generate. Furthermore, the SIMULATE program 
allows the user to select values of the original time series 
as starting values for the simulated series'. Output consists 


of the simulated series' and plots. 
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This completes the descriptions of the programs contained 
mn the Time Series Editor. The actual program listings are 
given in Appendix C and sample user sessions are shown in 
Appendix B. In the next chapter, an example is given of an 
entire time series analysis from plotting to diagnostic test- 


ing and forecasting. 
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770 EKANELE TIME SERIES ANALYSIS 


In this chapter a description is given of a complete 
analysis of a time series using the Box-Jenkins procedure 
and the Time Series Editor. The time series analyzed is 
series C (Chemical Process Temperature Readings) from Box 
and Jenkins [Ref. 2, p. 523]. This time series was selected 
because it is analyzed completely in Ref. 2 so that the 
user can compare the results found there with the results 
given by the Time Series Editor. 

The first step in the analysis of series C is to plot 
the time series. This is shown in Figure 3. The plot 
reveals rather wide fluctuations in the series but not the 
sort of explosive nonstationary behavior that would render 
a modeling attempt fruitless. The plot also reveals that 
the time series has a large amount of momentum (movements 
of the series tend to resist changes of ےت‎ This 
is characteristic of ARIMA models with one or two differences. 

Second, the autocorrelations, partial autocorrelations, 
mean, and variance of the series were estimated using AUTO. 
The numerical values are printed out in Table I and plots 
of the autocorrelations and partial autocorrelations are 
given in Figures 4 and 5. Figure 4 shows that the autocor- 
relations dampen out slowly in a near linear fashion. This 


15 ان‎ indication that the series is nonstationary and that 
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one or more differences are needed to make it stationary. 
The partial autocorrelation plot is not informative when 
the autocorrelations fail to dampen out rapidly. 

As suggested by the plots of original series and its 
autocorrelations, the program DIFF was executed to trans- 


form the series ا"‎ to the series tw, J where w = ለሃ, = 


በ. Yo: The values of the series tw) are tabulated in 


Maple Il. The autocorrelations and the partial autocorre- 
lations of the differenced series were then calculated and 
plotted. The values are given in Table III and the corre- 
lation plots are shown in Figures 6 and 7. The correlation 
plots suggest that the first differenced series is either 
ARMA (1,0) with the AR parameter near unity (the autocorre- 
lations of the first differences also dampen out slowly) 
or a second difference is required. Thus, two candidates 
are suggested: 

INNFA RIMA (1,1,0): (1-6,B) (1-B)y, = 0 +a 
and 2) ARIMA (0,270): (3-8) y, = ። 一 — 
For purposes of estimation, the second model was extended to 
include two moving average terms. SUG e EE 


often done to see if the estimated moving average parameters 


turn out to be near zero, thus confirming the tentative 


Identification. Thus, the two models entertained were: 
l) (1 9,B)Ay, SÉ a, ( ) 

2 2 
and 2) A vi 0 (1 9,B 9, ja, ( ) 


The next step is to calculate maximum likelihood esti- 


mates of the model parameters. The program ESTIMATE was 
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The autocorrelation and partial autocorrelations of the 
residuals were calculated and plotted to test the model. 


The correlations should appear to be estimates of a pure 


white noise process if the model is adequate. The model 
passed that test. A goodness-of-fit test was also performed 
to test the model. The chi-square value was = 2 8 ۲٣ 


19 df and a significance level of 0.3082. Thus, there is 
I trong evidence to suggest that the (1,1,0) model: 

(1-1 .8131B + 0.8131B8^)y, ar 11) 
is inadequate. 


Parameter estimates were also obtained for the ARIMA 


ü model using ESTIMATE. The parameter estimates were: 
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DN crore, the correlation plots of the residuals fail to 
suggest any inadequacy of the model. However, the chi- 


square lack of fit test yielded: 


3 = 28.74 with 18 df and a significance level of 0.0516. 
Thus, if the ARIMA(0,2,2) model were the correct model a 
chi-square value as large as 28.74 would occur by chance 
nma »robability of only 0.0516. There is some ground 
here for questioning the ARIMA(0,2,2) model. Because of its 
simplicity (parsimony is a very desirable feature for all 
models) and its better fit, the ARIMA(1,1,0) model, equation 
(11), was selected. 

The estimated model, equation (11), for series C was 
used to calculate forecasts and confidence limits for those 
forecasts. The forecasts were made for 13 periods into the 
future at origin t-20 using the observed values of the 
original time series as starting conditions. Figure 8 shows 
a plot of the forecasted values adjoined to the original 
series with EOS confidence limits about those forecasts. 

The forecasted values and the probability deviations are 
tabulated on page 59. The plot shows that the confidence 
limits are very wide for lags far into the future. 

Finally, program SIMULATE was utilized to generate two 
simulated series from eq. (11). Values of these two simula- 
tions are given in Table IV. Plots of those simulated 
Series are shown in Figure 9. The general shape of those 
mimos is like that of the original time series, thus con- 
firming the estimated model. 
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The keyboard printout of the user session that generated 
the analysis above is included in Appendix B. The entire 
session lasted approximately one and a half clock hours in- 
cluding checkout of plots and consumed less than one minute 


of CPU time. 
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vie SUMMARY AND RECOMMENDATIONS 


There is a growing need for an efficient unified collec- 
ESS of computer programs for aiding operations researchers, 
statisticians, economists, managers, and other people who 
must analyze time series. The Box-Jenkins methodology has 
Opened the doors of time series analysis to an expanding popu- 
lation of analysts. Many computational algorithms have been 
developed and are available in many forms. The problem has 
been the iterative nature of the Box-Jenkins' model building 
procedure. Such a procedure is straightforward but can be 
very time consuming. The Time Series Editor that has been 
described in this report provides a unified collection of 
computer programs in an interactive time-sharing mode that 
can aid the user in all phases of the model building procedure 
from the plotting of the time series to the forecasting of 
future values. The Editor does not develop new computational 
algorithms- Rather, it makes those that are available easier 
and faster to use. With its simple input requirements which 
SE. all prompted by written instructions, the Editor can 
easily be used by the most naive computer use. 

The Box-Jenkins methodology has been described in Chapter 
2 and descriptions of the modules of the Time Series Editor 
have been given in Chapter 3. These descriptions were given 


not as a substitute for study of the Box-Jenkins technique, 
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Eneas a communication device to explain to the potential 
user of the Time Series Editor what happens in the various 
stages of the computational process, and to explain what the 
Editor requests mean. An example time series analysis cover- 
ing all stages of the Box-Jenkins model building procedure 
was covered in Chapter 5. Appendices contain a User's Guide 
to the Time Series Editor, sample user sessions and program 
puc, and program listings. 

Although the Editor covers the entire model building pro- 
cess, as described by Box and Jenkins, from the plotting of 
the series to diagnostic testing and forecasting, much more 
could be added to improve the Editor's capability and utility. 
Listed below are several options that are recommended for 
addition to the Time Series Editor. 

jr. Extend the diagnostic ability to include a periodo- 

gram analysis or other tests related to the spectral 
analysis of time series. 

2. Include an option that will determine all roots of 

the characteristic equation and give the general 
solution to the autocorrelation function and to the 


eventual forecast function. 


3. Modify the FORECAST option to allow forecasts for 
seasonal nonstationary series'. 


4. Modify the ESTIMATE option to provide parameter esti- 
mates for seasonal stationary series'. 


Expand the univariate model building capability to 
linear transfer function model building. 


6. Expand the Time Series Editor to include multivariate 
models such as multiple regression. 
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APPENDIX A: USER SES GUIDE TO TIME SERIES EDITOR 


In order to use the Time Series Editor, the user must 
log onto CP, link to the disk storage area where the Time 
Series Editor resides, implement CMS, log into the general 
user and Time Series Editor disk areas, and enter the EXEC 
routine. It is also necessary to log out of the system at 
Ene completion of execution. The material which follows 
will enable the user to perform the above steps on the NPS 
CP/CMS system. Commands marked with an asterisk (*) are 
entered by the user (the asterisk itself is omitted). Those 
without an asterisk and those in all capital letters are 
system responses. Numbered sentences are comments which 
will not appear during an actual user session. The instruc- 
tions and system responses assume an IBM 2741 Input/Output 
Terminal. Some modifications may be necessary if other 
terminals are used. 


d Turn the terminal on, depress the RETURN key, and wait 
for the system to respond: 


CP-67 online xd.65 qsyosu 


ER: Depress the ATTN key. The roll bar will advance and 
cne keyboard will unlock. Then enter: 


login xxxxgnn 
3。 XXXX is the user's identification number, and nn is 
the terminal number (written on the right-hand-side 
of the terminal. For example, if the user's ID is 0655 
and the terminal number is 07 the input would be: 
Hogan 0655907 


4. The system will respond with the statement: 


ENTER PASSWORD: 
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The user then enters his password (most users at NPS 
have the password npg): 


*npg 
The system will then respond: 


ENER 4-DIGIT PROJECT NUMBER FOLLOWED BY 4-CHARACTER 
COST CENTER CODE: 


The user then enters: 
*aaaabbbb 


aaaa is the assigned project number and bbbb is the 
user's section designator or faculty code. 


The system will respond with the message of the day 
sueh as: 


CP/CMS HRS.. 0930-2200 (MON-THURS)..0930-1800(FRI) OUTPUT 
RETAINED 4 DAYS 

CEMS VERSION 3.2 

At this point the user is in CMS. He must then get in- 
EU. THIS is done by hitting the ATTN key. The 
system will respond: 


ER 


The user must then link to the TIME SERIES EDITOR. This 
is accomplished by entering: 


NE 19695 191 192 

The system will respond with the instruction: 
ENTER PASSWORD: 

The user then enters: 

ZI II 

The system will respond: 

SET TO READ ONLY 

The user must now implement CMS by entering: 
TDL cms 

The system will respond: 


CMS VERSION 3.2 
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Now the user must log into both his general user and 
the Time Series Editor area by entering: 


“login 191 


The system will respond with a message about the status 
of the file such as: 


ከከ በክ: (ንሪ FILES; 225 REC IN USE, 71 LEFT (OF 296), 
COS FULL (2 CYL) 
۳ 


The user should then enter the command: 


login 192 t,p 
R; 


The system will respond: 


929) R/O 
R; 


The user can then enter the Time Series Editor by 
issuing the command: 


*timeser exec 

The system will respond: 

MAA JE ENTERED THE TIME SERIES EDITOR 
ENE ASE RESPOND TO EACH INQUIRY 

etc. 

The user is then on his own, guided by the Exec routine. 
See the notes that appear at the end of this Appendix 
FOF additional information. Eventually, the user will 
be asked: 

DO YOU WANT TO TRY AGAIN? 

If a yes response is given another sequence will begin. 
A no will take the user out of the Time Series Editor. 


The system response will be: 


CONTROL RETURNED TO CMS 
R; 


The user then can log out by entering: 
Rep logout 


The system will respond: 
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AAN T Ee 001.03.37 TOTCPU= 001. 30. 1 4 
ከ ከ ከሚ 2171 1 7.19.083 ON 3/21/77 


p. .جس جح تھے‎ .፣ ስ ፡።፡።ክ túrn off his terminal and tear off 
lis output. 


A sample of the procedure is shown below: 
cp-67 online 20.65 11 


login 0655907 

ENTER PASSWORD: 

npg 

ENTER 4-DIGIT PROJECT NUMBER FOLLOWED BY 4-CHARACTER COST 
CENTER CODE: 

0986rk54 

CP/CMS HRS, MONDAY THRU THURSDAY,0930-2200,FRIDAY,0930-1800. 
OUTPUT RETAINED 5 DAYS 

READY AT 16.13.22 ON 03/22/77 

CMS VERSION 3.2 


offline read * 
P (191): 49 FILES; 225 REC IN USE, 71 LEFT (OF 296), 76% FULL 
O C YL) 


OFFLINE READ FILE FT02F001 
R; 

CP 

link 1969p 191 192 

ENTER PASSWORD: 

Kerr 

Soe nO READ ONLY 


ipl cms 
)ا‎ 135 VERSION 3.2 


191 یی 1 70 

P (191): 49 FILES; 225 REC IN USE, 71 LEFT (OF 296), 76% FULL 
(2 CYL) 

R; 

Megan 132 t,p 

T (192) R/O 

R; 

timeser exec 


UO HAVE ENTERED THE TIME SERIES EDITOR 


DO YOU WANT TO TRY AGAIN? 
n 
CONTROL RETURNED TO CMS 
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R3 


CB log 
SONNECT= 00.03.38 VIRTCPU- 000.00.34 TOTCPU= 000.01.30 
۶۰۰0 ٣ AT 13.48.56 ON 3/11/77 


NOTES: 


E BINS error In CP/CMS Cam be corrected by typing the 
character as many times as is required to backup and 
then type the correct values. For example, the command 
exec timesre could be corrected by typing two @ signs 
2۰۰۷۰۱۷ ٠ب‎ cne correct. letters er as follows: 


exec timesre@@er 


An entire line can be deleted by typing the character ¢. 


The user should exercise care before depressing the re- 
turn key. If an input is required and the return key 
is depressed before any character is entered, the user 


may get tossed out of the Editor. Before the return key 


is depressed, any errors can be corrected. After the 
return is hit there are only selected cases that can 
later be corrected. 


When asked to enter integer data, the user should always 


Meee justify the data within the allowed region. 


If at any time the user encounters a debug condition 
(probably caused by errors in entering data) and he 
cannot determine what needs to be done, he can get out 


of the debug condition by depressing the ATTN key twice. 


This puts the user back into CP from which he can link 
Dato the Time Series Editor or logout. 
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APPENDIX B 


SAMPLE USER SESSION AND PROGRAM OUTPUT 
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APPENDIX B: SAMPLE USER SESSION AND PROGRAM OUTPUT 


cp-67 online xd.65 qsyosu 


login 9655g07 
ENTER PASSWORD: 
npg 

ENTER 4-DIGIT PROJECT NUMBER FOLLOWED BY 4-CHARACTER COST CENTER CODE: 
0986rk54 

CP/CMS HRS, MONDAY THRU THURSDAY,0930-2200,FRIDAY,0930-1800. 
OUTPUT RETAINED 5 DAYS. 

FILES:- 01 RDR, NO PRT, NO PUN 

READY AT 16.13.22 ON 03/22/77 

CMS VERSION 3.2 


offline read * 
P (191): 49 FILES; 225 REC IN USE, 71 LEFT (OF 296), 76% FULL (2 CYL) 


OFFLINE READ FILE FTO2FOOI 
R; 

CP 

link 1969p 191 192 

ENTER PASSWORD: 

FETE 

SET TO READ ONLY 


ipl cms 
CMS VERSION 3.2 


login 191 
Ple do FILES; 225 REC IN USE, 71 LEFT (OF 296), 76% FULL (2 CYL) 


R; 

login 192 t,p 
T (192) R/O 
R; 

timeser exec 


YOU HAVE ENTERED THE TIME SERIES EDITOR. 

PLEASE RESPOND TO EACH QUERY WITH AN INPUT AT THE TERMINAL. 
ENTER ONLY THE FIRST LETTER FOR A WORD RESPONSE. 

ENTER NUMERICAL VALUES VIA FORTRAN FORMAT. 

TYPE INTEGER VALUES FOR NAMES STARTING WITH 1 THRU N. 

TYPE INTEGER VALUES FOR COUNTS. TYPE FLOATING 

VALUES WITH DECIMAL FOR ALL OTHER NAMES. 


DO YOU WANT A LIST OF OPTIONS? 
y 
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OPTION DESCRIPTION 


GENERATE 一 一 一 一 一 GENERATE ANY ARIMA TIME SERIES 

DTD — CALCULATE AUTOCORRELATIONS, PAUTOS, MEAN AND VARIANCE 

ID ------ PLOT A TIME SERIES 

۳۳ << سم اد‎ een TRANSFORM AND DIFFERENCE A TIME SERIES 

ላ ፣ 2---- CALCULATE MAX LIKELIHOOD ESTIMATES OF ARMA PARAMETERS 

851111. sees SIMULATE TIME SERIES FROM A GIVEN MODEL 

ESBECAST 一 一 一 一 一 FORECAST FUTURE VALUES, CONSTRUCT CONFIDENCE INTERVALS 


WOULD YOU LIKE MORE INFO? 
x 


ENTER OPTION YOU WANT INFO ABOUT. 
d 


+ THIS PROGRAM TRANSFORMS A GIVEN TIME SERIES WITH A LOG 

OR AN EXPONENTIAL TRANSFORMATION (OR NO TRANSFORMATION IF DESIRED), AND 
THEN TAKES SEASONAL AND/OR SIMPLE DIFFERENCES OF ANY SPECIFIED ORDERS. 
iT THEN OUTPUTS THE TRANSFORMED AND DIFFERENCED TIME SERIES. THIS 
PROGRAM IS USED TO ATTEMPT TO MAKE A SEASONAL OR A NONSTATIONARY TIME 
SERIES OF THE FORM THAT CAN BE HANDLED VIA BOX-JENKINS TECHNIQUES. IT 
ESESCEPRDIF EROM THE IMSL LIBRARY. 


DO YOU WANT TO TRY A SESSION? 
X 


ENTER THE LETTER FOR OPTION YOU WANT. 

p 

IS YOUR DATA IN FILE FT02F0012? 

y 

ENSSCUTION BEGINS... 

ENTER TITLE FOR PLOT. 

series c, original values. 

TINE SERIES PLOTS HAVE BEEN PRINTED OFFLINE. 
DO YOU WANT TO TRY AGAIN? 


Y 


ENTER LETTER FOR OPTION YOU WANT. 
a 
IS YOUR DATA IN FILE FTO2F001? 


Y 
SEXECUTION BEGINS... 


AUTOCORRELATIONS 


0.978 0.944 04902 0.854 0.802 0.748 0.692 0.635 0.570 3 
ASS 0215 00.359 105305 0.252. 0.201606150 0.098 0.047 -0.003 


PARTIAL AUTOCORRELATIONS 


0.978 -0.260 -0.157 -0.093 -0.058 -0.045 -0.012 -0.038 -0.022 -0.010 
-0.036 -0.041 -0.038 -0.024 -0.037 -0.027 -0.032 -0.070 -0.048 -0.024 
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MEAN= 22.9739 VARIANCE = 4.22273 


ENTER TITLE FOR PLOTS. 

series C, original values. 

YOUR AUTO AND PAUTO PLOTS HAVE BEEN PRINTED OFFLINE. 
PICK UP IN ROOM I140 UNDER YOUR USED ID NUMBER. 

DO YOU WANT TO TRY AGAIN? 


X 


ENTER LETTER FOR OPTION YOU WANT. 

d 

IS YOUR DATA IN FILE FTO2F001? 

4 

EXECUTION BEGINS... 

IS YOUR TIME SERIES SEASONAL? 

n 

ENTER NUMBER OF NONSEASONAL DIFFERENCES. 

| 

DO YOU WANT A LOG TRANSFORMATION? 

n 

DO YOU WANT AN EXPONENTIAL TRANSFORMATION? 

n 

YOUR TRANSFORMED TIME SERIES HAS BEEN STORED IN FILE FT03F001. 
DO YOU WANT A PRINTOUT OF THE FIRST TEN VALUES? 


Y 
THE FIRST 10 VALUES OF THE TRANSFORMED SERIES ARE: 
0.400009 07099991 ORO 0.0 0.0 
02199997 -0.100006 =O. 0995991 -0.300003 023799992 
DO YOU WANT TO PLOT AUTO ANDPAUTO OF TRANSFORMED DATA? 
Y 
*EXECUTION BEGINS... 
AUTOCORRELATIONS 


88:05 095.659 0.596 07442 0.380 0.318 0.262 0.186 0.139 0.144 
0.097 0.094 0.074 0.073 0.070 0.072 0.089 0.043 0.041 0 


PARTIAL AUTOCORRELATIONS 


۳٣۰٣٣٠۰۱0۰0710 -0.007 0.051 0.027 -0.019 -0.013 -0.080 0.020 0.117 
SIT 0-094 20.027 0.034 0.006 0.012 0.043 -0.121 0.060 -0.005 


MEAN=-.346667E-01 VARIANCE = 0.531982E-01 
ENTER TITLE FOR PLOTS. 
series C, lst order differenced data. 
YOUR AUTO AND PAUTO PLOTS HAVE BEEN PRINTED OFFLINE. 
PICK UP IN ROOM 1140 UNDER YOUR USED ID NUMBER. 


DO YOU WANT TO TRY AGAIN? 
y 


ENTER LETTER FOR OPTION YOU WANT. 
d 
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IS YOUR DATA IN FILE FTO2F001? 

y 

EXECUTION BEGINS... 

IS YOUR TIME SERIES SEASONAL? 

n 

ENTER NUMBER OF NONSEASONAL DIFFERENCES. 

2 

DO YOU WANT A LOG TRANSFORMATION? 

n 

DO YOU WANT AN EXPONENTIAL TRANSFORMATION? 

11 

YOUR TRANSFORMED TIME SERIES HAS BEEN STORED IN FILE FTO3FOOl. 
DO YOU WANT A PRINTOUT OF THE FIRST TEN VALUES? 


D 
THE FIRST lO VALUES OF THE TRANSFORMED SERIES ARE: 

20.300018 20-1 0.0 00 202199997 

8209999] 0.000015 =O. 200012 207099991 0.199982 
DO YOU WANT TO PLOT AUTO AND PAUTO OF TRANSFORMED DATA? 
y 
BIECUTION BEGINS... 
AUTOCORRELATIONS 


-0.079 -0.065 -0.022 -0.063 0.013 -0.018 0.049 -0.052 -0.124 0.122 
ከከከ ፒፒ በህ. -0.077 0.029 -0.011 -0.058 0.171 -0.101 -0.013 -0.020 


PARTIAL AUTOCORELLATIONS 


-0.079 -0.072 -0.135 -0.094 -0.022 -0.051 0.022 -0.060 -0.145 0.094 
-0.143 0.023 -0.089 -0.005 -0.039 -0.074 0.130 -0.105 -0.014 -0.053 


== .26/866E-02 VARIANCE = 0.198143E-01 
ENTER TITLE FOR PLOTS. 
series C, 2nd order differenced data. 
YOUR AUTO AND PAUTO PLOTS HAVE BEEN PRINTED OFFLINE. 
PICK UP IN ROOM I140 UNDER YOUR USER ID NUMBER. 


DO YOU WANT TO TRY AGAIN? 

E 

ENTER LETTER FOR OPTION YOU WANT. 
e 

*IS YOUR DATA IN FILE FTO2FO001? 
y 

EXECUTION BEGINS... 

ENTER NUMBER OF AR PARAMETERS. 

0 

ENTER NUMBER OF MA PARAMETERS. 

2 

ENTER NUMBER OF DIFFERENCES. 

2 

* LENGTH OF TIME SERIES = 226 

O AR PARAMETERS 2 MA PARAMETERS 2 DIFFERENCES. 
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ESTIMATED AR PARAMETERS (UNDIFFERENCED FORM) ARE: 
PHI (1) = 2 900 


BET (2) = -1.0000 


ESTIMATED MA PARAMETERS ARE: 
۳۱۲۶۰ (1) = 0.1382 


THETA (2) = 0.1300 
MA CONSTANT= 0.0 WHITE NOISE VARIANCE= 0.189515E-01 


DO YOU WANT TO PLOT AUTO AND PAUTO OF RESIDUALS? 
ሃ 
EXECUTION BEGINS... 

AUTOCORRELATIONS 


00.034 -0.134 -0.093 -0.012 -0.043 0.004 -0.063 -0.125 0.088 
ES 0.065 -0,098 0.025 0.004 -0.055 0.157 -0.089 0.006 -0.006 


PARTIAL AUTOCORRELATIONS 


Bebe -0.090 0.001 -0.055 -0.019 -0.071 -0.143 0.087 
IO -0,104 -0,013 -07025 -0.090 0,123 -0.123 -0.021 -0.025 


MEAN=- .124975E-02 VARIANCE - 0.189499E-OI 


ENTER TITLE FOR PLOTS. 

series C, residuals of (0,2,2). 

YOUR AUTO AND PAUTO PLOTS HAVE BEEN PRINTED OFFLINE. 
PICK UP IN ROOM 1140 UNDER YOUR USER ID NUMBER. 

DO YOU WANT TO CHECK RESIDUALS? 

y 

LXECUPION BEGINS... 

ENTER IP, NO OF AR PARAMETERS (DIFFERENCED FORM). 

0 

ENTER IQ, NO OF MA PARAMETERS. 

2 

* CHI SQUARE LACK OF FIT VALUE = 28.74 DF = 18 SIGNIFICANCE = 0.0516 


DO YOU WANT TO TRY AGAIN? 

y 

ENTER LETTER FOR OPTION YOU WANT. 
e 

TSE YOUR DATA IN FILE FTO2FOOL? 
ሃ 

EXECUTION BEGINS... 

ENTER NUMBER OF AR PARAMETERS. 
1 

ENTER NUMBER OF MA PARAMETERS. 
0 

ENTER NUMBER OF DIFFERENCES. 

1 


37 


AIDA ur (TIT AR 











Ww a 1 NËT 57 ROT 
و ار‎ SIT OT 


és TUA ٠ 
26 7ہ‎ 7 SA o. 


. 19 D AS AM 2 


PNT OM 8 


LENGTH OF TIME SERIES = 226 
] AR PARAMETERS O MA PARAMETERS tore PE RENCES. 


ESTIMATED AR PARAMETERS (UNDIFFERENCED FORM) ARE: 


PHI(1) = 221 
PHI(2) = -0.8131 
MA CONSTANT= 0.0 WHITE NOISE VARIANCE- 0.178396E-Ol 


DO YOU WANT TO PLOT AUTO AND PAUTO OF RESIDUALS? 
Di 
EXECUTION BEGINS... 

AUTOCORRELATIONS 


0.010 0.009 -0.053 -0.009 0.058 0.021 0.076 -0.020 -0.089 0.130 
20:095 0,083 -0,056 0,042 0.006 -0.041 0.169 -0.080 -0.001 -0.010 


PARTIAL AUTOCORRELATIONS 


ES. 0.009 -0.054 -0.008 0.060 0.018 0.074 -0.016 -0.088 0.140 
=0.104 0.069 -0.048 0.041 0.003 -0.030 0.157 -0.087 0.013 -0.025 


MEAN=-.897935E-02 VARIANCE = 0.177590E-01 


ENTER TITLE FOR PLOTS. 

series C, residuals of (1,1,0). 

YOUR AUTO AND PAUTO PLOTS HAVE BEEN PRINTED OFFLINE. 
PICK UP IN ROOM I140 UNDER YOUR USER ID NUMBER. 

DO YOU WANT TO CHECK RESIDUALS? 

ሃ 

EXECUTION BEGINS...» 

ENTER IP, NO OF AR PARAMETERS (DIFFERENCED FORM). 

1 

ENTER IQ, NO OF MA PARAMETERS. 

0 

* CHI SQUARE LACK OF FIT VALUE = 21.53 DF = 19 SIGNIFICANCE = 0.3082 


DO YOU WANT TO TRY AGAIN? 

Y 

ENTER LETTER FOR OPTION YOU WANT. 

f 

IS YOUR DATA IN FILE FTO2F001? 

M 

EXECUTION BEGINS... 

ARE YOUR PARAMETERS IN DIFFERENCED FORM? 
Y 

ENTER NO. AR PARAMETERS (DIFF. FORM). 
1 

ENTER AR PARAMETER PHI(1). 

058131 

AR PARAMETERS ARE: 

0.8131 


28 


2225596 
224.2594 


0.966487 
191838 


3.80429 
4.80092 


ARE THESE OK? 

X 

ENTER NUMBER OF MA PARAMETERS. 

O 

ENTER OVERALL MA CONSTANT 

O 

ENTER THE NUMBER OF DIFFERENCES. 

1 

ENTER INDEX FOR FORECAST ORIGIN VIA 13. 
020 

ENTER FORECAST LEAD TIME VIA 12. 

13 

ENTER SIGNIFICANCE LEVEL, ALPHA. 

2 

MA CONSTANT = 0.0 NUMBER OF DIFFERENCES = 1 

FORECAST LEAD TIME = 13 ALPHA = 0.200 

ARE THESE OK? 

y 

AR PARAMETERS IN UNDIFFERENCED FORM ARE: 
18131 

-0.8131 


THE FORECASTS ARE: 


23.1561 227 22.7964 22.6653 
2284719 22.4014 22.3440 2232974 
2.2286 2222035 22.1831 

80.0% CONFIDENCE LIMITS FOR FORECASTS ARE: 

0.172546 03547271 0.556689 36153 
ን ን ነሪ 1236521 5035 1.74057 
2.08959 2425639 2.41304 

PSI WEIGHTS FOR FORECAST ERROR ARE: 

0 +0 2.47423 3.01150 3.44889 
4.09326 4.32822 4.51927 4.67461 
4.90362 4.98 712 5. 405501 


WHITE NOISE VARIANCE = 0.18127E-01 
DO YOU WANT TO PLOT FORECASTS? 


y 
*DO YOU REALLY WANT TO PLOT FORECASTS? 


y 
YOUR TIME SERIES AND FORECASTS ARE PLOTTED OFFLINE. 


DO YOU WANT TO TRY AGAIN? 
Y 


ENTER LETIER FOR OPTION YOU WANT. 

S 

IS YOUR DATA IN FILE FTO2F001? 

y 

*EXECUTION BEGINS... 

ENTER NUMBER OF AR PARAMETERS (UNDIFFERENCED FORM). 
2 
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18.781552 
10.639331 
20256069507 
20.241806 


ENTER ESTIMATED AR PARAMETER PHI(1). 


00 


ENTER ESTIMATED AR PARAMETER PHI (2). 


-0.8131 


۶۰۰1 


AR PARAMETERS ARE 
ARE THESE OK? 


ሃ 


ENTER NUMBER OF MA PARAMETERS 


0 


ENTER OVERALL MA CONSTANT 


0 


ENTER ESTIMATED WHITE NOISE VAR 


020178 


DO YOU WANT TO INPUT STARTING VALUES? 


n 
START (1) = 
START (2) = 


STARTING VALUE (S) OK? 


ሃ 


ENTER RANDOM NUMBER SEED (BETWEEN O AND 1). 


07659376 


ENTER NUMBER OF VALUES YOU WANT TO SIMULATE VIA 13. 


856260951 
19.261063 
20.347610 
2095297957 


065 


DO YOU WANT MORE THAN ONE SIMULATED SERIES? 


18.631083 
., 1 
20166321 
20.174332 


M 


HOW MANY SERIES DO YOU WANT (LE.9)? 


2 


DO YOU WANT A PRINTOUT OF THE FIRST 20 VALUES? 


ሃ 


SIMULATED VALUES ARE: 
18.640320 
16.963882 
9821411 
202216507 


DO YOU WANT TO PLOT SIMULATED VALUES? 


y 


EXECUTION BEGINS... 
BETTER TITLE FOR PLOT. 


simulations with model (1,1,0), series C. 
TIME SERIES PLOTS HAVE BEEN PRINTED OFFLINE. 
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DO YOU WANT TO TRY AGAIN? 


n 


CONTROL RETURNED TO CMS 


R; 
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Residuals of (0,2,2) model of series C. 
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